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Abstract 

This paper studies the leading chiral corrections to heavy quark-diquark symmetry predictions for 
doubly heavy baryon semileptonic decay form factors. We derive the couplings of heavy diquarks to 
weak currents in the limit of heavy quark-diquark symmetry, and construct the chiral Lagrangian 
for doubly heavy baryons coupled to weak currents. We calculate chiral corrections to doubly 
heavy baryon zero-recoil semileptonic decay for both unquenched and partially quenched QCD. 
This theory is used to derive chiral extrapolation formulae for measurements of the doubly heavy 
baryon zero-recoil semileptonic decay form factors in lattice QCD simulations. 
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I. INTRODUCTION 



Motivation for this work comes from the SELEX experiment's recent observation of states 



which have been tentatively interpreted as doubly charm baryons 



2|, l3| . While the masses 



and the hyperfine splittings of the observed states are consistent with theoretical predictions 



from some quark models 



6|, lZI, many other aspects are 



J] and quenched lattice QCD 
difficult to understand theoretically. Therefore model-independent predictions for doubly 
heavy baryon properties would be helpful. Chiral and heavy quark symmetries are useful 
approximate symmetries of QCD for predicting low-energy properties of hadrons containing 
a single heavy quark. Heavy quark-diquark symmetry relates heavy mesons to doubly heavy 



antibaryons. The 



splitting 



5 



13 



leavy quark-diquark symmetry prediction for the doubly charm hyperfine 



wrote an effective La- 

n 



14( 1 is within 25% ~ 30% of the quenched lattice QCD calculation. More 
observables are needed to see if heavy quark-diquark symmetry can be applied to charm and 
bottom. In this paper we will study the leading chiral corrections to heavy quark-diquark 
symmetry predictions for zero-recoil doubly heavy baryon semileptonic decay form factors. 

Heavy quark-diquark symmetry relates mesons with a single heavy quark (Qq) to an- 
tibaryons with two heavy antiquarks {QQq). Savage and Wise 
grangian for heavy quarks and diquarks using heavy quark effective theory (HQET) [9] and 
obtained a prediction for doubly heavy baryon hyperfine splittings. HQET considers only 
the scales Aqcd and ttlq and is formulated as an expansion in Kqcd/^q- Hadrons with 
two or more heavy quarks, such as quarkonia or doubly heavy baryons, are characterized 
by the additional scales rnqv and mQv"^, where v is the typical velocity of the heavy quark 
within the bound state. The correct theory of hadrons with two or more heavy quarks is 



Non-Relativistic Quantum Chromodynamics (NRQCD) 



121. Refs. 



ij] derived 



2 



effective Lagrangians for heavy diquarks within the framework of NRQCD and corrected a 
factor of 2 in the prediction for hyperfine sphttings of Ref. 8|]. In this paper we will use 
these techniques to derive the couplings of heavy diquarks to weak currents. 

For studying the low energy dynamics of hadrons containing diquarks, it is useful to 



build an effective theory which incorporates t 



15 



16 



le relevant symmetries of QCD. Heavy hadron 



17| has heavy hadrons, Goldstone bosons, and 



chiral perturbation theory (HHxPT) 
photons as its elementary degrees of freedom and incorporates chiral and heavy quark sym- 



metries. Ref. 



181] derived an effective Lagrangian for doubly heavy baryons incorporating 



heavy quark-diquark symmetry and used this theory to study low energy strong and elec- 
tromagnetic interactions of doubly heavy baryons. In this paper, we will apply this theory 
to calculate chiral corrections to semileptonic decay at zero-recoil. 

Lattice gauge theory is a useful tool to address the nonperturbative dynamics of quarks 
and gluons. First principles numerical calculations of doubly heavy baryons from lattice 
QCD can be compared with experiment to test our understanding of QCD. Lattice cal- 
culation of decay matrix elements of heavy mesons and doubly heavy baryons can help 
determine the reliability of heavy quark-diquark symmetry for charm and bottom hadrons. 
HPQCD and MILC collaborations are planning to calculate semileptonic decay form factors 
of heavy mesons in the near future. In lattice QCD, they do simulations with unphysical 
sea quark masses which need to be extrapolated to the physical values. Chiral Lagrangian 
in quenched and partially quenched theory is useful for the chiral extrapolations of lattice 
data. In this paper, we will extend the chiral Lagrangian with heavy quark-diquark symme- 
try to include effects of quenching and partially quenching, and use it to derive formulae for 
the chiral extrapolation of zero-recoil semileptonic decay of double heavy baryons in lattice 
QCD simulations. 
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The rest of this paper is organized as follows. In Sec. [TTl we use NRQCD to derive the 
couplings of heavy diquarks to weak currents consistent with heavy quark-diquark symmetry. 
In Sec. Illlt we construct the chiral Lagrangian for doubly heavy baryons coupled to weak 
currents and calculate the doubly heavy baryon zero-recoil semileptonic decay matrix ele- 
ments. We calculate the chiral corrections to doubly heavy baryon zero-recoil semileptonic 
decay using xPT. In heavy quark-diquark limit the corrections vanish. In Sec. llVt we extend 
the chiral Lagrangian with heavy quark-diquark symmetry to the partially quenched theory 
and derive formulae for chiral extrapolation of doubly heavy baryon zero-recoil semileptonic 
decay form factors in lattice QCD simulations. We compare chiral corrections to the zero- 
recoil doubly heavy baryon semileptonic decay form factor in PQxPT with the corresponding 
corrections in xPT. In Sec. |Vl a brief summary is given. Some useful formulae are collected 
in the Appendix. 



II. COUPLINGS OF HEAVY DIQUARK TO WEAK CURRENTS FROM VN- 
RQCD 

NRQCD is the nonrelativistic effective theory which describes the dynamics of non- 
relativistic heavy quarks. In NRQCD there are four important energy scales, uiq, the 
typical momentum of the heavy quarks within the bound state mqv, the typical kinetic 
energy of the heavy quarks rnqv'^ and Aqcd- The total momentum of the heavy quark 
field is taken to be the sum of the label momentum p and the residual momentum k as 
Ptotai = p(~ ^qv) + k(~ rnQv"^). vNRQCD is an effective theory for NRQCD which has 
a consistent v expansion constructed by Luke, Manohar and Rothstein 12|. The leading 
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order vNRQCD Lagrangian to 0{v^) for the QQ sector for b, c antiquarks is 

_ i V V - Y^tf ^y^y^^ -V - y^tf ^y'^y"^ f^y" + ■ ■ ■ (1) 

f=b,c p,q p,q 

where the elhpsis represents high order terms. Xp is a nonrelativistic antiquark field with 
flavor / which annihilates an antiquark. The kinetic energy Dq and momentum p of the 
heavy quarks are 0{mQV^) and 0{mQv), respectively. is SU{3) color generators for 3 
representation, B is chromomagnetic field, Dq and D are the time and spatial components 
of the gauge covariant derivative. In order to derive an effective Lagrangian for diquarks we 



follow the methods of Ref. 



14| . We use the spin and color Fierz identities 

1 . . . . 1 



SasSp-y = --{o-'e)ap{ea')^s + -eape^s ^ (2) 

m (mn) 

to project the potential onto the 3 and 6 channels and to decompose the antiquark bilinears 
such as XpX-p iiifo operators with spin and spin 1. In Eq.([3]) The Greek letters denote spin 
indices, the a* denote the Pauli matrices and e = zcr^. In Eq.(l3]), the Roman letters denote 
the color indices and the matrices d^J"^^ are symmetric matrices in color space defined by 



n 



5™ 5" m = n 



For bb and cc, the diquark state must be in either(3)c(3)5 or (6)^(1)5 by the Pauh principle. 
However, there is no restriction on the color and spin for mixed flavor. After a Fourier 



5 



transform to obtain position space potentials, the Lagrangian can be rewritten as 

/=f),c V q / V P 

+•••, (5) 

where the eUipsis reflexes to terms in which the diquark is in the 6 of color, V^^\r) = — |^ 
and V^^\r) = The first line includes the gauge boson fields, the kinetic terms for the 
antiquarks and the leading spin symmetry breaking interactions which generate hyperfine 
splittings. The next three lines are the quartic terms of diquark fields. In order to get 
the interaction for doubly heavy baryon semileptonic decay, we need the Lagrangian for 
diquarks of different flavor. The diquark fields for cc and bb can be introduced using the 



Hubbard-Strantonovich transformation as in Ref. 



141 ]. We add to the Lagrangian: 



f=b,c V q J 



■^-q^J 



Tr-J2e^--'-^{x':,he^{4)i], (6) 



2 

p 



where m = 1,2,3 is the color index. The = e*P''"^(x{p)A:ecr(x^)z is the 
spin-1 // field, t" = Ep e*P''"^(X-p)fceo-(Xp)« is the spin-1 be field, and the T'™ = 
e«p-r£r|M.(^T^c ^^^gTj^^fe -j^ -g ^YiQ spin-0 be field. The quartic terms of antiquark fields cancel 
and we are left with the following action 

C = 

f=b,c \ p q 

\ p q 

V p q / 

The Feynman rules for the couplings of the diquark fields to the pairs of antiquarks from 

the interaction terms in Eq.(I7]) are shown in Fig. [H Integrating out the fields T^, T™ and 
T^™ will give us back the original NRQCD Lagrangian. Integrating out the antiquark fields 
Xn and X-n "^ill yield an effective action for the diquark fields T^, and T^"^. The be 



Lp "'^^^ A— p 

fields kinetic terms are the same as those of cc fields in Ref. 



14| . with the exception that 



the diquark reduced mass is now 



^ /^Q = • 8) 

z mbrric 



The coupling of the diquarks to weak currents is obtained by evaluating graphs in Fig. [21 
The fiavor changing weak currents are given by Jq^d = C7^(l — 75)6. By matching onto 
NRQCD, the weak current is rewritten by JjJ-^ = x'^K^'^^ '~ 5'"cr*)x'' + ■ • • where the ellipsis 
reflexes to higher order terms. For the lowest order weak currents from Fig. [2] we find the 
couplings of diquark to weak currents. The Lagrangian for the weak interactions of diquarks 



7 



J- J. , i J. 



-P, 7, 



FIG. 1: Feynman rules for the coupling of the diquarks T™, and T^"^ to antiquarks. The Greek 
letters denote spin indices and the Roman letters refer to the color indices. 



T.r 




FIG. 2: One loop diagram contributing to the coupling of composite anti-diquark fields and weak 
currents. 



IS 



(9) 



where Jw is the weak current. In the last line we expand the diquark fields to the lowest 
order and write the lagrangian in terms of local currents. The ■ ■ ■ denotes the excited diquark 
states contribution and the factor t] can be interpreted as the spatial wavefunction overlap 
of the ground state initial diquark with the ground state final diquark system, which is not 



predicted by symmetry, i.e. tj = J (i^r0^^(r)0cc(i') e**' '". Our results differ from Ref. 



;2lj by a 



different relative sign of the coupling to the zeroth component of the weak current and that 
of the spatial component of the weak current. 



III. DOUBLY HEAVY BARYON SEMILEPTONIC DECAY 

In the heavy quark limit, the lowest mass states of the spin-0 heavy meson, P, and the 
spin-1 heavy meson, P*, are degenerate, and therefore combined into a single field. Ha, 

Ha = P:-a + Pa, (10) 

where a is an SU{3) flavor anti-fundamental index and the a are the Pauli matrices. Simi- 
larly, the ground state degenerate doubly heavy baryon doublet. Tip ^, consists of a spin-1/2 
doubly heavy baryon, S^^^, and a spin-3/2 doubly heavy baryon, S* 



Ta,ip = V2 ( S* + -J=^a,^ (tI,p ] , (11) 



1 „ 

where S*^^^ and Ea^.y are the spin-| and spin-| fields, respectively. The i = 1, 2, 3 are diquark 
spin indices and /3 = 1, 2 are light antiquark spin indices. The field S* .^^ obeys the constraint 
'^*a,ii3 '^^7 — 0- Heavy quark-di quark symmetry relates the heavy mesons to doubly heavy 
antibaryons. The effective Lagrangian for doubly heavy baryon and heavy meson ground 



18| 



state doublets in the heavy hadron rest frame was constructed in Ref. 

£ = TT[ni{tDo)bani>] - g TTiniHk a ■ + ^Tr [7^1 Ha cr'] , (12) 

where 'Ha,^p = Ha^afi + ^a,i/3 is a 5 X 2 matrix field which transforms as tensor products of 
the five component field (which corresponding to the two heavy quark spin states and 
the three diquark spin states) and a two-component light antiquark spinor. Here the index 



^ In this section, the T is the doubly heavy baryon field, although T is diquark field in the previous section. 



/i takes on values between 1 and 5, a, (3 = 1 or 2, and i = 3,4, or 5. The first term is the 
kinetic term, the second term is the coupling to axial current and these two terms respect 
the SU{5) heavy quark-diquark symmetry. The third term is the hyperfine splitting term 



which is the leading heavy quark symmetry breaking. The S^,^ are the 5x5 matrices 

(13) 



y Tjk J 



k 



and (X^)jk = —i^ijk- is the hyperfine splitting of the charmed mesons which is related 
to the hyperfine fplitting of the doubly charm baryon by heavy quark-diquark symmetry. 
Aba is the spatial part of the axial vector current. 

The relevant terms of Eq. f|T2l) for our calculation are 

+ T'+iDoT'^ - g'T';T'p,ap,p ■ A. (14) 
The composite fields are decomposed as 

(Ty^ = ^ , (15) 

where S*^ and Sq, are the spin-| and spin-| fields, respectively. Here, we have suppressed 
the fiavor indices. Tj/3 and Tj/j are the ground state doublet of ccq and bcq , respectively, and 
the diquarks are in spin-1 channel. is the spin-| ground state of bcq where the diquark 
is in spin-0 channel. The hyperfine splitting of the doubly heavy baryons are related to 
that of the charmed mesons by heavy quark-diquark symmetry, m=* — = jAh and 



m=*^ — ms^^ = (j^j^^ where the /iq is the reduced mass of diquark be. 
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Eq{T3] describes the low energy strong and electromagnetic interactions in the doubly 
heavy baryon rest frame in which the doubly heavy baryon four-velocity is conserved (up to 
0{Aqcd/^q) corrections). For a process such as the weak decay, in which the initial and 
final baryons have different four-velocities, the covariant representation of baryon fields is 
needed. However, for studying the zero-recoil semileptonic decay, in which the doubly heavy 
baryon four-velocity is conserved, it is possible to work in the baryon rest frame. Therefore 
we are allowed to match the weak currents coupling of diquark onto a current operator in 
the effective theory for doubly heavy baryons. The Lagrangian for semileptonic decays of 
doubly heavy baryon is 



The couplings of the doubly heavy baryon to the weak currents is obtained by demanding 
that its diquark index on the doubly heavy baryons couple to the weak currents as in Eq.Q. 

Expanding out Eq. fll6l) in terms of doubly heavy baryon fields, we find the following weak 
current matrix elements in terms of the initial and final spinor. 



(16) 



^bcq) — V "^q' 



( 




(17) 




(18) 




(19) 




(20) 






(22) 
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FIG. 3: One-loop contributions to the doubly heavy baryon semileptonic decay. 

where {ua, u'^, Uip) and {ua,Uij3,u'^) are nonrelativistic spinor for initial and final state, re- 
spectively. These form factors are predicted by heavy quark-diquark symmetry. The tree 



2Q| 



level zero-recoil doubly heavy baryon semileptonic decay matrix elements agree with Ref. 
up to an overall minus sign^. At tree level, the 6i are zero and we confirm that the 6i van- 
ish at any order in the heavy quark limit. They will receive corrections from heavy quark 
symmetry breaking effects. Obtaining those symmetry breaking corrections that come from 
chiral loop corrections in Fig. [3] is the goal of this paper {6[ and 62 are not calculated in this 
paper) . 

The pion-baryon vertex in Fig. [3] is generated from the pion-baryon interaction terms in 
Eq.f ll4l) and the weak current vertex is from the expanded weak currents terms in Eq. f|T6l) . 
There are also loops for the wavefunction renormalization which are not shown here. The 
chiral loop corrections to semileptonic decay form factors are given in the Appendix and a 
table of 6i with different values of fi is given in Table [I] for both q = u,d and q = s. The 
PQxPT corrections will be discussed in the next section. Here we work in the isospin limit, 
so ruu = rrid. We take the the reduced mass of diquark be to be fiq = ~ rUc. We 



^ They give results for Ebc while we give the results for S^g. To obtain their result to the lowest order, we 
should use the decomposition of the ground state doublet as Tip — \/2 — '^'^}3.y'^a,'yj and then we 
agree up-to an overall sign. 
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oi[u, a) 


62{u, d) 
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04(^U, u ) 




6q{u, d) 
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-0 05 


-0 05 
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-0 01 


-0 09 




di{s) 


02(,Sj 










= 500 MeV 


-0.01 


-0.01 


-0.14 


-0.14 


-0.16 


-0.16 


;U = 1500 MeV 


-0.0002 


-0.001 


-0.13 


-0.12 


-0.13 


-0.13 



TABLE I: Chiral correstions to doubly heavy baryon containing au,d or s semileptonic decay form 
factors for /u = 500 MeV and /i = 1500 MeV. 

n 

choose g = 0.6 [23] and we set f-K = fk = fri = 130 MeV. The experiment value for Ah is 
140 MeV. In the numerical calculations we use the values of Goldstone masses of ~ 140 
MeV and ^ 500 MeV and = y — 2_ — *l from the . Then only one parameters can 
be varied, the renormalization scale /x. In the chiral perturbation theory calculations the 
logarithmic /i-dependence from loops is cancelled by /i-dependent counterterms we have not 
included here. We vary /i from 500 MeV to 1500 MeV to obtain an estimate of uncertainty 
due to the unknown counterterm contributions. We find that, for the doubly heavy baryon 
containing antiup or antidown quark, the first and second form factors get very big chiral 
corrections ~ 25% — 32% and the fifth and the sixth get small correction ~ 5%. The form 
factors for doubly heavy baryon containing antistrange quark obtain negative corrections 
and the first and second form factors are least sensitive to chiral corrections ~ 1%. It would 
be interesting to see if the observed deviations from heavy quark-diquark symmetry in either 
experiments or lattice results agree with the chiral corrections. 
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IV. PQxPT RESULTS 



The EFT techniques with heavy quark-diquark symmetry for doubly heavy baryons can 
also be applied to doubly heavy systems simulated on the lattice. It is useful to extend the 
xPT results to include lattice artifacts such as quenching and partial quenching. Here we 
will focus on partially quenched chiral corrections to the doubly heavy baryon zero-recoil 
semileptonic decay form factors, which provide the formulae needed for chiral extrapolations 
in lattice calculations. The major modification to the calculation in the previous section is 
that there is a modified propagator of the Goldstone mesons. 

In a lattice calculation, the sea quark masses are often different from valence quark masses. 
In partially quenched QCD the sea quark masses are different from the valence quark masses, 
while in quenched QCD, the sea quark contributions are absent. To reproduce the lattice 
unphysical quarks in field theory, fictitious ghost quarks are added to cancle the valence 
quark contributions and then we are left with the effects of sea quarks. This is equivalent 
to the partial quenching artifact in lattice QCD simulations which use different masses for 
sea and valence quarks. In the limit of rrisea = f^vaience, QCD is recovered. PQ^PT is useful 



for lattice extrapolations to physical quark masses. The PQy 



meson sector is described by the Lagrangian 



24 



25 



26 



21, 



28 



'T anc 



29 



30, 



QyPT pseudoscalar 



3l|, 



f2 7 

C = y str (S^S+a^S) + - str (m^S^ + mjs) + a<i,9^<l>o9^$o - 



(23) 



where the field E is defined by 



exp 



(24) 
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and the meson fields appear in the f/(6|3) matrix, 



M 



(25) 



The operation str() in Eq. fl23l) is a supertrace over flavor indices, i.e., str(y4) = €aAaa, 
where = 1 for a = 1 — 6 and ea = — 1 for a = 7 — 9. The M and M matrices contain 
bosonic mesons, while the x ^-^id matrices contain fermionic mesons (one ghost quark 



with one sea/ valence quark). The quark mass matrix is deflned as 



22| 



nig = diag{mu, m^, m^, rrij, mi, rrir, niu, nid, nis). 



(26) 



We work in the isospin limit for both the valence and sea sectors, so m„ = rrid and mi = 
mj. We take the strange sea quark to be same as the valence quark, = m^. From 
the lowest order PQxPT Lagrangian, the meson with quark content qq' has the mass of 



m 



qq' f 



-4-(mg + mg/). The PQxPT propagators of the off-diagonal mesons have t 



le usual 



Klein-Gordon form. The flavor neutral propagator can be conveniently written as \ 



22 



Qah^ — ^a^abPa + 'Pab {Pa, Pb, Px) i 



(27) 



where 



Pb 



Tab {A,B,C) 



{mla - m]j) {ml 



T - ^bb + "^^ 
ml^) 



P 



X 



''aa 



^Ib) 



?2 



A + 



[m 



bb 



q^ — + 
- m]^ {ml^ 



mz 



m 



X) 



^Ib 



-B 



{m\ 



m^j- ) [m 



ml 



-c 



m 



bb) 



where = | (m^^- + 2m: 



(28) 



In terms of the 5x2 fleld, the partially quenched Lagrangian for doubly heavy baryon 
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and heavy meson is 22j 



+ ^ (T^ts ■n(j)+a (n^HM) + a' (n^n) str(A<), (29) 

where M. is the mass operator defined by = | {C,fnq^ + ^^^g^^) Therefore the baryon 
mass sphttings in PQ^PT are given by 

^ccaa' = rn^ , — rri':^— = m~* — m-* = —aimai — nia) , 

^'^11 ^ccq> "ccq "ccq' ^ccq \ 1 9/ ' 

^hrnn' = t^i" , — 1^",. = 1^"* — t^i"* = —cr(m„> — m„) , 

'"^11 ^bca' ^bcq ^h^ijl ^1.7,7, \ Q <?/ ' 



-bcq' -^bcq -j^^ 

3 



Kcqq' = "^H- - m^,,q = -^H - (y{m^, - m,j , 



where we take the the reduced mass of diquark he to be iiq = ^ nic, a is the couphng 

constant in the mass operator in the Lagrangian in Eq. (|29l) . We choose g = 0.6 23| and / = 
130 MeV. The experiment value for Ah is 140 MeV. With the values of the SU{3) splitting 
of the ground state D mesons — Tno = —a {iris — ^u) ~ 100 MeV, pion mass value 



m. 



^ = jiirriu + md) and kaon mass value mi = jj{ms + m^) , we obtain —a = ^^_^2 ji- For 
the mass of r] we use the SU{?>) prediction, = \J "''^+4^ ^^ Then only three parameters 
can be varied, the valence pion mass = "^mu, the sea pion mass m-^sea = "^jj, and 

the renormalization scale /z. In the chiral perturbation theory calculations the logarithmic 
/i-dependence from loops is cancelled by /i-dependent counterterms we have not included 
here. We vary /i from 500 MeV to 1500 MeV to obtain an estimate of uncertainty due to the 
unknown counterterm contributions. In the heavy quark limit, the logarithm /x dependence 
goes away and we obtain the tree level predictions in previous section. Calculating the 
PQxPT loop diagrams, such as Fig. [3l we find the corrections for the six form factors which 
we list in Appendix. 
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200 300 400 500 600 700 200 300 400 500 600 700 

m^sea(MeV) m^sea(MeV) 



FIG. 4: A5n {q — u) as a function of Tn-j^sea foi^ different values of fn-j^uai- The width of the bands 
is the results of varying ^ between 500 MeV and 1500 MeV. 

In the Figs. H] and [5] we plot the relevant corrections to the form factors A5j = 5iPQ — 6iy. 
in terms of m^^ga for different values of m^^a^ for g = m and q = s. For each value of 
iTT'nvai = 140 MeV, 280 MeV and 420 MeV, we let the m,rva/ range from mT^vai up to the mass 
of eta-strange, m^^ = rriss ~ 700 MeV. The bands correspond to varying /i from 500 MeV 
to 1500 MeV, which is chosen to be the same for both xPT and PQxPT. From the plots, it 
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FIG. 5: A6n {q = s) as a function of m-jrsea but independent of different values of rriT^yai- The width 
of the bands is the results of varying fi between 500 MeV and 1500 MeV. 

is easy to see that the partially quenched chiral corrections reproduce the chiral corrections 
when the sea quark mass goes to the physical valence quark mass. The corrections of 
the doubly heavy baryon containing a strange quark do not depend on the valence quark 
masses. As demonstrated by Figs. H] and [5l the partially quenched chiral non-analytic 
corrections relevant to the chiral corrections are very insensitive to /x- dependence. The 
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relevant corrections to 5i, ^2, and ^4 are affected most by sea quark masses ~ 10% — 40% 
and increase with increasing rusea values. The ^5 and 6q are very insensitive to partial 
quenching effects, receiving ~ 1% — 8% small relevant corrections from partial quenching, 
and are insensitive to sea quark masses, while the other four are not. It would be interesting 
to test the effects on lattice. Quark-diquark symmetry predictions should be more robust. 
Deviation from quark-diquark symmetry must be attributed to short distance effects. 

V. SUMMARY 

In this paper we have used NRQCD to derive the couplings of heavy diquarks to weak 
currents with heavy quark-diquark symmetry. We constructed the chiral Lagrangian for 
doubly heavy baryons coupled to weak currents and calculated the tree level predictions 
for doubly heavy baryon semileptonic weak decay form factors. We calculated the chiral 
corrections in both unquenched and partially quenched theory. The partially quenched loop 
corrections to semileptonic decay form factors are given in the Appendix. The formulae will 
be useful for chiral extrapolation of doubly heavy baryon zero-recoil semileptonic decay form 
factors in lattice QCD simulations. The chiral corrections of the six form factors range from 
— 15% to 30%. S5 and Sq for doubly heavy baryon containing antiup or antidown quark receive 
small corrctions ~ 5% and Si and S2 for doubly heavy baryon containing antistrage quark 
receive smaller correction ~ 1%. The partially quenched chiral non-analytic corrections 
relevant to the chiral corrections are very insensitive to //-dependence. The ^5 and 5q receive 
small relevant corrections ~ 1% — 8% from partial quenching, while the other four are not. 
It is interesting to see if the size of the corrections in experiments and lattice calculations 
match the chiral derivations of the calculations. 
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APPENDIX: xPT AND PQxPT FORM FACTORS 

Here we give the detailed one-loop xPT and PQxPT corrections to zero-recoil semilep- 
tonic decay. In the calculation we use MS scheme and have not included the countert- 
erms. In the chiral corrections to the six form factors / is pion decay constant, is the 
mass of the Goldstone boson in the one-loop diagram and C^^ is a factor which comes 
from SU (3) Clebsch-Gordan coefficients in the couplings. For loops with charged pions we 
have = = 1, for loops with neutral pions Ci° — C22 = |, for loops with kaons 
= Cf = 1 (i =1 or 2), and for loops with 77 mesons C^^ = C^2 = | and C^g = |. 

i,q' ^ ' 

4 20 

- ^H^lcqq', ^ccqq', m, /^) " ^HKcqq'^ Kcqq'^^i^ /^)) > (32) 

,2 



^3xiQ) = J2 ^ll' (47r/)2 ^~27^^^'"'''^'' -Kcqq'^rUi, /x) + —I{Abcqq', A^cgg', TTli, /x) 



'± 20 

- ^H^lcqq', -Kcqq'^ /^) " 1^1 {^Icqq' ^ ^ccqq' , /x)) , (33) 



<^4x(g) = J2 ^m'J^^^-^§f^(~^bcqq'^ ^'^l' ^ ^) + ^^-^bcqq': Kcqq': ^i, fl) 

4 20 

- ^I{^bcqq', Accqq',mi, fj,) - — /(Afocgg', A*^^^,, TUi, //)) , (34) 
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g"^ f A 5 \ 

<^5x(9) = Cig, ^^^^^^ \^--I{-Algg,, -A^^^^,, 7^^, z^) ■ -I{Ahcgg', ^00^^'^^^, /x) j ,(35) 



4 11 

- ^li^bcgq', -Al^gg,,mi, ii) - —I{Abcgq', Accgg' , rrii, ii)) . (36) 

The formulae for xPT corrections to doubly heavy baryon semileptonic decay form factors 
are given in terms of the following functions, 

7(0,0, m,/x) = 0, 

/(A, A, m, li) = -6A2 In + 4 (m^ - SA^) + 8mAF + 4 (A^ - m^) F' , 

7(Ai, A2, m, 11) = -2 (A? + A1A2 + A^) In + 4 (m' - A? - Ai A2 - A^) 



where 

^(,)^^'-v^(f— -'(7*P)).W<1 . (38) 
y/x^ — l\n {x + V^^^ — 1) , kl ^ 1 
and F\x) is the first derivative of x. 

Next we give the same corrections for the partially quenched case. 

g"^ 1 8 

5iPq{(1) = -(J^i-gH^bcgg': Accgg', m,,/, //) - -I{Al^g,, A*^^,, rUgg, , //)) 

g'=j,l,r ^ ■' ' 

g"^ 1 8 
^ (47r/)2 (~9^(^''^gg' ^'^^99' "^99' '"^ ~ 9^(^6cgq> A*cqq, m^g, m^g, n)) ,(39) 

^2PQ(g) = /Z^N2 (^-^(^fc^gg^ ^C'^^g'' - ^H^bcgq', A*^^^,, m^g/, /x) 

- ^H^bcqqh Accqq', rUgqi , fj) - — /(A^^^^,, A*^^^,, mgg/, /x)) 

1 4 ^ 

~l~ "(477^)2^ ^27^^ ^55' ^qi?) A'') ~ ^-^(Afecgg) ^ccqqi '^qqi ^qqi f^) 

4 20 

- 27^(^6cgg> Accqq, ^gg, m^g, //) - — /s:(A^,gg, A*,^^, rugg, rugg, //)) , (40) 
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q'=j,l,r ^ ^> 

4 20 

+ (4^7)2 (~ ^-^(^^^99' "^99' ^99' /^) + ^^{Abcqq, ^ccqq, rUgg, TUgg, //) 

4 20 
- —K{Al^gg, -A*^gg, 7^^,, 7^^^, z^) " — ^(A^,^^, AccQg, m^^, m^^, //)) , (41) 

^^8 5 

^4PQ(g) = ^ /^^f\2i~'^H-^bcqq'^ ^ccqq',mgg>, fJ,) + ^/(-A^^^^,, A*^^^,, mgg/, /x) 

4 20 

-^H^bcqq', ^ccqq', ^ggi , /J,) - — /(Afocgg', A*^^^,, Ulgg, , n)) 
(^2 g g 

_ ±j^f. . v_20 

2Y^\^bcqqi ^ccqqilT^qqilTT'qqi fJ-) 



K{Abcqq, Accqq, TUgg, TUgg, fJ.) - —K{Abcgg, A^^g, m^g, rUgg, H)) , (42) 



f A 5 \ 

kpQ{q) = J] Ti^ff)^ ( ~9^(~^6cg?'> -^ccqqM ^9?', /^) " 0^(^6059', ^ccqq' , ^gq', /^) ] 

g"^ f A 5 \ 

+ ^^^^y [^-^^i-^lcqq^ ' Kcqq^ ^qq, ^qq, /^) " gH^bcqq, ^ccqq, TUgg, m^g, n) j^^) 

g2 ^ ^ 

(^6PQ(g) = Yl /4^^p (~^-^(~^&cgq^ -Kcqq'^^qq'^ /^) " ^^-^bcqq': ^ccqq' : TUgg,, fl) 

4 ^ 11 

- -^H^bcqq', -Kcqq'^ '>^qq' ^ /^) " -^H^bcqq', ^ccqq', rUgg' , fJ,)) 

g2 Q 4 

+ (47j-y^2 (~'2J-^(~'^bcqqJ ~^ccqqJ ^qqi ^qqi ^) ~ l^^{~^hcqqi ^ccqqi ^qqi ^qqi ^) 

4 11 

- ^K{Ahcgq, -A*^gg, TUgg, rrigg, fx) - —K{Abcgg, A^gg, TTlgg, Ulgg, fl)) , (44) 

where the function K{Ai, A2, m, m, //.) which arises from the hairpins is given by 

K{Ai, A2, rria, rrib, n) = Tab (H^i, A2, rua, n)., /(Ai, A2, m^, /x), /(Ai, A2, mx, A*)) • (45) 
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